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ON  THE  SIGNIFICANCE  OF  A WALL  EFFECT  IN 
ENCLOSURES  WITH  GROWING  FIRES 

Leonard  Y.  Cooper 
Abstract 

This  paper  studies  the  significance  of  a wall  effect 
that  has  been  observed  during  the  growth  stage  of  enclo- 
sure fires.  Relative  to  the  two-layer  phenomenon  which 
tends  to  develop  during  such  fires,  the  effect  has  to  do 
with  the  near-wall  downward  injection  of  hot  upper  layer 
gases  into  the  relatively  cool  uncontaminated  lower 
layer.  It  is  conjectured  that  these  observed  wall  flows 
are  buoyancy  driven,  and  that  they  develop  because  of  the 
relatively  cool  temperatures  of  the  upper  wall  whose 
surfaces  are  in  contact  with  the  hot  upper  layer  gases. 
For  a growing  fire  (growth  proportional  to  tm;  t being 
time  and  m 0)  in  an  enclosed  compartment,  an  analysis 
of  the  conjectured  mechanism  for  the  wall  flow  leads  to  a 
time-dependent  solution  for  the  ratio  of  wall  flow  mass 
ejection  rate  from  the  upper  layer,  m^,  to  the  fire  plume 
mass  injection  rate  _to^  the  upper  layer,  m^.  The  solution 

indicates  that  in  practical  fire  scenarios  m /m  can  be 

w p 

of  the  order  of  "several  tenths"  even  prior  to  the  time 
that  the  upper  layer  interface  has  dropped  to  an  eleva- 
tion midway  between  the  ceiling  and  fire.  In  other 
words,  the  results  of  the  analysis  indicate  the  impor- 
tance of  taking  the  wall  effect  into  account  in  two-layer 
zonal  analyses  of  enclosure  fire  phenomena. 
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1 . INTRODUCTION 


The  zone  approach  to  mathematical  modeling  has  received  much  attention  in 
the  analysis  of  enclosure  fire  environments.  In  recent  years  several  such 
zone  models  have  been  developed  to  predict  environments  in  rooms  of  fire 
origin  as  well  as  in  single  and  even  multiple  adjacent  spaces.  Because  of  the 
fact  that  they  have  each  tended  to  be  developed  for  specific  purposes  and 
classes  of  problems,  there  is  quite  a variation  from  model  to  model  in  the 
amount  of  attention  and  detail  paid  to  the  simulation  of  different  physical 
phenomena  that  come  into  play.  In  this  regard,  however,  even  the  most  sophis- 
ticated of  the  models  have  tacitly  ignored  a commonly  observed  and  potentially 
important  wall  effect  which  can  lead  to  significant  features  of  enclosure  fire 
environments  not  heretofore  simulated. 

After  describing  the  wall  effect  in  question  this  paper  will  then 
identify  the  conditions  under  which  the  effect  is  significant.  Finally,  the 
paper  will  propose  a general  purpose  algorithm  for  including  the  wall  effect 
in  an  overall  zonal  model  analysis. 

2.  A DESCRIPTION  OF  THE  WALL  EFFECT  AND  A MEASURE  OF  ITS  SIGNIFICANCE 

The  wall  effect  under  consideration  is  first  described  and  analyzed 
within  the  context  of  a generic  fire  scenario  which  develops  within  an 
enclosure  lacking  any  significant  ventilation  openings  (any  leakage  from  the 
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enclosure  is  assumed  to  occur  near  the  floor) . Later,  the  relation  of  the 
wall  effect  to  fire  environments  developing  in  enclosures  with  open  doors  or 
windows  will  be  addressed. 

A sketch  of  the  significant  features  of  all  the  conjectured  wall  effects 
to  be  described  below  is  presented  in  figure  1. 

A fire  starts  a distance  H below  the  ceiling  of  an  enclosure  of  floor 
area,  A,  and  releases  energy  at  the  rate  Q(t),  where  t is  time  from  the  effec- 
tive instant  of  ignition.  Because  of  their  elevated  temperature,  the  products 
of  combustion  from  the  fire  are  driven  upwards  by  buoyancy  forces,  and  a 
turbulent  plume  is  generated.  All  along  the  axis  of  the  plume,  quiescent 
ambient  air  is  entrained  laterally  into  and  mixed  with  the  plume  gases.  The 
plume  eventually  impinges  on  the  ceiling  and  spreads  radially  outward  along 
this  surface  forming  a relatively  thin,  turbulent  ceiling  jet.  The  ceiling 
jet  gases  redistribute  themselves  across  the  entire  ceiling  area  of  the  enclo- 
sure, start  to  fill  the  upper  enclosure  volume,  and  eventually  submerge  all 
the  continuous  ceiling  jet  flow  activity.  Below  this  jet  activity  the  now 
reduced  momentum  plume  gases  form  a relatively  quiescent,  elevated 
temperature,  upper  gas  layer  which  continues  to  increase  in  depth  as  the  plume 
gas  upward  filling  process  continues  in  time.  The  elevation  above  the  fire  of 
the  lower  interface  of  this  gas  layer  is  designated  by  Z^(t). 

Except  for  regions  near  the  bounding  ceiling  and  wall  surfaces,  it  is 
reasonable  to  characterize  the  upper  layer  as  being  uniform  in  composition 
with  an  absolute  temperature,  Tu(t)  and  with  concentrations  of  products  of 
combustion  i,  C^(t)  (units  of  product  per  unit  mass  of  bulk  upper  layer 
mixture) . 
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The  ceiling  surface  temperatures  are  generally  different  from  Tu. 
surface  is  heated  from  the  initial  ambient  temperature,  T , by  convection  from 
the  high  temperature  ceiling  jet  gases  and  by  radiation  from  the  fire's 
combustion  zone.  Later  into  the  fire,  radiation  exchanges  from  the  bulk  upper 
layer  gases  at  Tu  and  from  other  bounding  surfaces  of  the  enclosure  may  also 
become  significant. 

In  general  the  temperatures  of  the  wall  surfaces  are  also  different  from 
Ty.  The  radiant  heating  mechanism  of  these  surfaces  is  similar  to  that  for 
the  ceiling.  However,  except  for  possible  locations  of  strong  ceiling  jet- 
wall  surface  flow  interactions  (i.e.,  where  wall  surfaces  are  close  enough  to 
the  plume-ceiling  impingement  point) , the  vigorous  convective  heat  transfer  at 
the  ceiling  is  generally  absent  at  the  walls.  Thus,  one  can  anticipate  that 
the  wall  surfaces  will  be  heated  from  the  initial  temperature,  T , at  a signi- 
ficantly lower  rate  than  the  ceiling  surface. 

It  is  conjectured  and  seems  to  have  been  observed1 » that  the  condition 
of  relatively  cool  walls,  at  absolute  temperature,  T , bounding  the  elevated 
temperature  upper  layer  gases,  at  Tu,  would  lead  to  the  development  of  a 
classic,  two-dimensional,  downward-directed,  natural  convection  or  buoyancy- 
driven  boundary  layer  flow  along  the  upper  portion  of  the  vertical  wall 
surfaces.  The  phenomenon  would  occur  in  this  basic  form  away  from  locations 
of  vigorous  ceiling  jet-wall  interactions.  The  boundary  layer  flow  would 
originate  near  the  ceiling,  increase  in  mass  and  momentum  flux  with  decreasing 
elevation,  and,  at  the  upper  layer  - lower  layer  interface,  be  injected  into 
the  lower  layer  of  the  enclosure  as  a downward-directed  wall  jet  of  upper 
layer  gases.  The  injected  gases  would  have  a bulk  temperature  lower  than  the 
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upper-layer  gases  but  higher  than  the  relatively  cool  and  uncontaminated 
lower-layer  gases. 

Once  in  the  lower  layer,  the  downward-directed  wall  jet  of  the  now 
upwardly-buoyant,  product  of  combustion  laden  gas  would  be  buoyed  back  upward 
and  away  from  the  wall  to  either  mix  with  and  contaminate  the  lower  layer  air 
or  to  entrain  additional  (i.e.,  in  addition  to  the  fire  plume)  lower  layer  air 
into  the  upper  layer. 

Whatever  the  eventual  disposition  of  the  wall  flow  contaminents,  the 
basic  wall  flow  phenomenon  could  significantly  alter  the  rate  of  development 
of  life-threatening  conditions  in  the  enclosure.  If,  for  example,  the  flow  is 
buoyed  upward  as  a wall  plume  and  re-enters  the  upper  layer,  then,  having 
entrained  additional  lower  layer  air,  the  depth  of  the  upper  layer  will  grow 
more  rapidly  than  it  would  have  otherwise,  albeit  at  a reduced  temperature  and 
product  concentration.  On  the  other  hand,  if  the  wall  flow  is  mixed  with  the 
lower  layer  gas,  then  the  net  effect  would  be  a contamination  of  the  lower 
layer,  less  rapid  growth  in  upper  layer  depth  and  at  increased  temperature  and 
concentration. 

When  the  mass  flux  leaving  the  upper  layer  in  the  wall  flow  is  large 
enough  to  significantly  retard  the  velocity  of  descent,  dZ^/dt,  of  the  upper 
layer  - lower  layer  interface  or  to  lead  to  significant  increases  in  tempera- 
ture and/or  product  of  combustion  contamination  of  the  lower  layer,  then  this 
flow  must  be  taken  account  of  in  a mathematical  simulation  of  the  overall 
enclosure  environment. 
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Let  nip  be  the  mass  flux  of  the  plume  into  the  upper  layer,  i.e.,  at  Z^. 
Also,  let  m^  be  the  total  mass  flux  of  the  wall  boundary  layer  flow  ejected  at 
from  the  upper  layer.  Then,  consistent  with  the  above,  an  analysis  of  the 

enclosure  environment  which  does  not  include  an  accounting  of  the  wall  flows 

• • • • 

can  lead  to  reliable  results  only  when  m « m . Thus  when  m /m  <<  1 the 

w p w p 

wall  effect  can  be  neglected,  and  when  this  is  not  the  case  it  must  be  taken 
into  account.  This  will  be  referred  to  below  as  the  Wall  Effect  Criterion. 

It  is  noteworthy  that  the  described  wall  flow  phenomenon  would  not  neces- 
sarily occur  only  in  rooms  of  fire  origin.  Indeed,  it  can  also  play  an 
important  role  in  the  redistribution  of  the  products  of  combustion  in  spaces 
which  are  adjacent  to,  or  communicate  with  rooms-of-f ire-origin.  Here,  the 
fire  generated  plume  in  the  above  room-of-f ire-origin  scenario  would  be 
replaced,  for  example,  by  an  analogous  upper  doorway  plume  which  was  generated 
by  inflowing  combustion  products.  "Smoke  layering”  in  the  adjacent  space 
would  be  initiated  in  the  usual  way.  Also,  at  least  until  relatively  late 
into  the  fire  when  radiative  heating  of  adjacent  space  wall  surfaces  could 
become  substantial,  these  surfaces,  being  only  weakly  heated  by  the  conjec- 
tured buoyancy  driven  wall  flows,  would  likely  remain  close  to  their  original 
ambient  temperature. 

The  above  discussion  is  from  the  perspective  of  the  early  growth  stage  of 
the  fire  when  a definitive  (Tu  - Tw)  temperature  differential  exists,  when  no 
significant  differential  between  the  lower  wall  temperature  and  the  lower  gas 
temperature,  T^,  (both  of  which  are  initially  at  Tj  has  developed,  and  when 
the  overall  wall  flow  activity,  as  described  earlier,  is  conjectured  to  occur. 
As  time  moves  on,  however,  a substantial  difference  between  TL  and  the 
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(probably)  relatively  greater  lower  wall  temperature  (heated  primarily  by 
radiation)  may  begin  to  develop.  At  least  near  the  floor  of  the  enclosure 
this  would  (probably)  tend  to  lead  to  an  upward-directed,  lower  wall  flow.  It 
is  also  possible  that  such  a lower  wall  flow  would  become  relatively  strong  at 
a time  when  the  downward-directed,  upper  wall  flow  starts  to  decrease  in 
strength  by  virtue  of  a reduction  in  (T  - Tw) . Under  such  a circumstance, 
the  nature  of  the  interaction  of  such  combined,  counter-current  wall  flows 
could  be  important  in  an  adequate  description  of  the  overall  enclosure 
environment;  the  early  growth  stage  description  of  a dominant  upper  wall  flow 
would  have  to  be  revised. 


An  analysis  of  the  combined  wall  flow  effect,  which  may  come  into  play  at 
later  stages  of  the  fire,  is  outside  the  intended  scope  of  the  present  paper. 
The  combined  flow  effect  will  not  be  discussed  further  except  to  note  that  a 
study  has  been  recently  initiated  on  the  significance  of  the  phenomenon  in 

o 

ventilated  enclosure  fire  scenarios  during  steady  state  conditions. 

3.  SOME  PROPERTIES  OF  THE  WALL  FLOW 

Refering  to  the  wall  flow  depicted  in  figure  1,  the  total  mass  flux, 
m(x) , at  a distance  x = H - Z below  the  ceiling  (Z  above  the  fire)  can  be 
estimated,  for  Prandtl  Number  (Pr)  = 0.72  and  according  to  whether  the 
boundary  layer  is  laminar  or  mostly  turbulent,  from 


m(x) 

P 


1.70  y Gr 


1/4 


2/5 

0.102  y Gv 
00  x 


Laminar:Grx  < 0.5(10^) 
Turbulent  :Grx  > 1.5(10^) 


(1) 
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where  P is  the  total  length  of  the  perimeter  of  the  enclosure,  Grx  is  the 
Grashoff  number 


Gr 


x 


T v2 
00  00 


(2) 


and  where  dimensionless  values  of  Tu,  Tw,  and  Z have  been  defined  as 


T /T  ; 
u 00 


<t>  = T /T  ; 

Tw  w 00 


Z/H 


(3) 


The  above  estimates  for  m(x)  were  established  from  results  presented  in 
references  4-6.  The  Grx  bounds  for  a laminar  or  mostly  turbulent  boundary 
layer  are  from  reference  7.  The  Bousinesq  approximation  was  used  in  develop- 
ing eq.  (1),  and  will  be  used  throughout  the  analysis  to  follow.  All  gas 
properties  will  be  taken  as  those  of  air  at  standard  ambient  conditions.  The 
wall  temperature  is  taken  to  be  uniform. 

From  the  origin  of  the  boundary  layer  to  some  x above  the  interface  (0  < 
x < H - Z^),  the  average  rate  of  heat  transfer  per  unit  area,  q" , to  the  wall 
can  be  estimated  from 


q” 


(4) 


where 


Average  Nusselt  Number  = 

1/4  5 

0.48  Gr  , Laminar 
x ’ 

2/5 

0.0184  Gr  , Turbulent 
x ’ 


6 


(5) 
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and  where  k is  the  thermal  conductivity  of  the  gas  mixture.  Below  the  inter- 
face the  heat  transfer  rate  will  be  less  than  q" . 


The  mixing  cup  temperature,  T , of  the  wall  flow  at  position  x is  defined 
by 


T 

u 


T 

m 


q"?x 

mC 

P 


where 


is  the  specific  heat. 


(6) 


Using  Pr  = 0.72,  the  results  of  eqs.  (1),  (4)  and  (5)  in  eq.  (6)  lead  to 
the  following  estimate  for  Tm: 

T - T {0.39,  Laminar 

-T— = ) (7) 

u w ( 0.25,  Turbulent 

This  last  result  indicates  that  Tm  is  independent  of  x and  somewhat  closer  to 
Tu  than  to  Tw.  In  particular,  eq.  (7)  is  valid  at  the  interface  position 
x = xi  = H - Zi< 


The  kinematic  momentum  flux,  K,  of  the  wall  flow 


K = /ou2dy  (8) 

where  u is  the  downward  velocity  and  y is  the  distance  from  the  wall  surface, 
can  be  estimated  from  results  of  integral  boundary  layer  analyses  in  refer- 
ences 6 and  7.  Thus,  for  Pr  = 0.72 
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(9) 


K 


2.25v2  Gr3/4/x, 
00  x 


Laminar 


0.0366V2  Gr9/10/x,  Turbulent 
00  x 

Finally,  the  order  of  magnitude  of  the  downward  velocity  will  be  of 
interest  in  later  discussion.  From  results  in  references  4 and  6 it  is 
possible  to  obtain  the  result  that  for  Pr  = 0.72  the  maximum  velocity  in  the 
boundary  layer,  umax,  for  both  laminar  and  mostly  turbulent  flow  can  be  esti- 
mated from 


u 

max 


0.55V 

00 


(10) 


4.  IMPLICATIONS  OF  THE  WALL  EFFECT  CRITERION 


The  mass  flux  of  the  plume  at  the  upper  layer  - lower  layer  interface  can 

O Q 

be  estimated  from  * 


m 


0.21p„g1/2H5/2(l  - 


X )l/3Q*l/3;5/3  !/3 

r ' xo  i T 


(11) 


where 


- VH,  * 


Q(t)/Q  . Q. 


p»cpto//2h5/2 


(12) 


and  where  Q0  is  a characteristic  energy  release  rate,  and  (1  - Xf)  is  the 
fraction  of  Q which  effectively  acts  to  heat  the  plume  gases  and  ultimately 
drive  the  plume's  upward  momentum.  X is  approximately  the  fraction  of  Q lost 
by  radiation  from  the  combustion  zone  and  plume.  For  hazardous  flaming 

O 

fires,  Xf  is  typically  in  the  range  0.3  - 0.4.  In  real  fires,  the  simple 
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description  of  Eq.  (11)  has  been  shown  to  yield  excellent  results  above  the 
combustion  zone  when  the  "point"  location  is  taken  to  be  at  the  lowest  eleva- 
tion of  the  combustion  zone  or  plume  above  it  where  "free  unrestricted", 
lateral  entrainment  of  ambient  air  is  possible.® 


In  order  to  invoke  the  Wall  Effect  Criterion,  eqs.  (1)  and  (11)  are  used 

to  estimate  m /m  with  the  result 
w p 


m /m 
w p 


oGr  (P/H) 


0.2l(l-Xr)1/3Q*1/3^/3(gH 


3,  2 
/v„ 


U/2  .1/3 

J 


(13) 


where 


Gr  = Gr  (x 
x 


x.  = H-Z  ) = (<j)  -l)(l-6)(l-£  )3gH3/ 
i i u i 00 


(14) 


6 = U - l)/(<j>  - 1) 
w u 


(15) 


a = 1.70. 


_ 1 / /. 


a = 0.102,  n = 2/5  if  Gr  > 1.5(1010),  Turbulent 


(16) 


To  apply  eq . (13)  to  the  dynamic  fire  scenario  discussed  in  the  last  section 

it  is  now  necessary  to  estimate  the  values  of  t.(t),  d>  (t)  and  <t>  (t). 

i u w 

In  the  absence  of  any  wall  effect,  and  from  ignition  (^  = 1)  up  to  the 
time  that  the  interface  elevation  starts  to  approach  the  elevation  of  the  fire 
(i.e.,  becomes  small),  the  solutions  for  C^(t)  and  4>u(t)  have  been 
previously  obtained^  for  enclosure  fires  where  Q(t)  ~ tm  (arbitrary  m 0) , 
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and  where  any  pressure  relieving  leakage  from  the  enclosure  is  primarily  near 
the  lower  portions  of  the  enclosure  boundary.  Under  such  circumstances 


♦u-  { l-(m  + 3)eT3(,n+1)/(m+3)/[3(m+l)(l-C1)]r1 


(17) 


where 


e - (l-Xc)l[(n,+3)/3]mQ*}  2/('t3)/(l-»r)(’tl)/(,t3)  (18) 

t = 3[(l-Xr)0*]1/3(tH3/2g1/2/A)<m+3)/3/(m+3)  (19) 

and  where  Qq  is  explicitly  defined  as 

Qq  - (Q/tm)(H3/2g1/2/A)"m  (20) 


Also,  Xc  is  defined  as  the  fraction  of  Q which  is  instantaneously  transfered 

to  the  internal  bounding  surfaces  of  the  enclosure.  XCQ  would  include  all 

heat  transfer  to  surfaces  by  both  convection  and  radiation,  and  Xc  is 

g 

typically  in  the  range  0.6  - 0.9.  For  magnitudes  of  energy  release  rates  and 

1 2 

enclosure  sizes  of  practical  interest  , e <<  1 and  can  be  estimated  from 

lim  5 - (1  + 0.140t)“3/2[1  + 0(e)]  (21) 

e-K) 

for  5 and  m which  include  the  ranges  0.5  1 and  0 _<  m 2 . (For  a 

more  definite  idea  of  the  utility  of  the  eq.  (21)  approximation  refer  to 
figure  2 of  reference  12.)  Unless  noted  otherwise  all  further  estimates  in 
this  paper  will  only  be  applicable  when  the  above  approximation  yields  an 
acceptable  estimate  for  £^(t). 
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Using  eqs.  (17)  and  (21),  and  eliminating  the  explicit  dependence  of 
on  time  leads  to 


lim  (<f>  -1) 
e+0  U 


e(m+3) 

3(m+l)(l-ci) 


3(m+l)/  (m+3) 

+0(e)2 


(22) 


Using  eq.  (22)  in  eq.  (14)  and  rewriting  eq . (13)  eventually  leads  to 


aGrnX, 


m /m 
w p 


0.21xJ/2^/3 


m_ 

(m  + 

_ 


3) 


(23) 


Gr  = Gr(^i;  X^  m,  6) 


(l-6)(m+3)  f )2  ( + ' 1 ) 

3 ( m+ 1 ) VA  \ 0.140  / 


3(m+l) 

(m+3) 


(24) 


where  a and  n are  given  in  eq.  (16),  and  where  X^  and  X2  are  defined  according 
to 


(25) 


4.1  Applying  Results  When  Wall  Temperatures  Are  Close  to  Ambient 


During  the  early  growth  stage  of  the  fire  the  increase  of  Tw  above  its 
initial  value,  T^,  will  often  develop  at  a much  slower  rate  than  will  the 
increase  of  T . Thus,  for  the  purpose  of  evaluating  the  significance  of  the 
wall  effect,  it  is  reasonable  to  investigate  the  implications  of  eqs.  (23)  and 
(24)  under  the  condition 


6 « 1 


(26) 
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According  to  these  equations,  and  for  in  = 0,  1.0,  and  2.0,  figure  2 presents 


plots  of  C . as  a function  of  Xi  with  (m  /m  )/X_  as  a parameter,  under  the 

1 W p L 

condition  6=0.  Thus,  for  0 m < 2.0,  and  for  a specified  X^  fire  scenario, 

the  plots  of  figure  2 can  be  used  to  determine  the  value  of  below  which  the 

wall  effect  starts  to  become  significant  (i.e.,  m /m  becomes  larger  than, 

w p 

say,  several  tenths). 


The  results  of  figure  2 will  be  used  to  evaluate  the  significance  of  the 
wall  effect  during  test  runs  of  two  full-scale  enclosure  fire  test  programs 
reported  in  references  1 and  2.  As  mentioned  earlier,  the  wall  effect  seems 
to  have  played  a role  in  those  experiments. 

The  first  evaluation  refers  to  a constant  16.2  kW  enclosure  fire  reported 
in  reference  2.  The  fire  was  an  acetylene  diffusion  flame  generated  from  a 
gas  burner  source  2.06  m below  the  ceiling  and  in  the  center  of  a square, 
"fully"  enclosed  space  (designed  to  leak  from  below)  with  sides  of  length 
3.64  m.  Tests  on  the  burner  indicated  a 50%  radiation  loss  from  the  sooty 
flame  (X^  = 0.5).  Total  loss  to  the  ceiling  and  walls  was  estimated  to  be 
84%  (Ac  = 0.84).  In  terms  of  previous  definitions,  and  using  the  values 
T^  = 29 4°K , = 1.2  kg/m3,  = 1.5  (10~5)  m2/s,  Cp  = 240  cal/(kg°K),  and 

g = 9.8  m/s  here  and  in  all  later  calculations,  the  test  run  is  found  to  be 
characterized  by 


m = 0 , e = 0 .0036 , X^ 


1.4(109),  X2  = 4.0 
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• • 

Now  assume  the  wall  effect  to  be  significant  when  m /m  > 0.3  , i.e., 

w p 

(m  /m  )/X„  > 0.075  . Then  figure  2 for  m=0  leads  to  the  conclusion  that  the 
w p z 

wall  effect  must  be  taken  account  of  once  drops  to  0.62,  i.e.,  once  the 
upper  layer  thickness  exceeds  0.78  m. 

The  remaining  evaluations  will  refer  to  test  runs  reported  in  reference  1 
in  a 2.36  m high,  fully  enclosed  two  to  three  room  space  (designed  to  leak 

O 

from  below),  where  one  of  the  rooms,  having  an  area  of  14.0  m and  designated 
as  the  burn  room,  contained  a centrally  located  methane  burner  2.12  m below 
the  ceiling.  The  one  to  two  rooms  adjacent  to  the  burn  room  were  corridor- 
like in  configuration,  and  the  total  multi-room  test  space  area  ranged  from 

O 

40.6  - 89.6  m . For  test  runs  with  all  rooms  freely  connected  (i.e.,  full 
open  doorways),  data  indicated  that  after  an  initial  time  interval  a single 
room,  two  layer  mod  had  some  applicability  in  analyzing  the  dynamic  conditions 
within  the  overall  space.  It  is  therefore  reasonable  to  utilize  the  present 
room-of-f ire-origin  results  in  an  evaluation  of  the  significance  of  the  wall 
effect  under  consideration.  As  mentioned  earlier,  the  existence  of  the  wall 
effect  was  observed  during  the  course  of  both  the  reference  1 and  reference  2 
test  programs. 

The  first  evaluation  of  the  above  test  series  will  involve  the  smallest 

O 

constant  fire,  25  kW,  and  the  largest  test  space  area,  89.6  m of  the  test 
program.  This  configuration  had  a total  wall  perimeter  of  approximately 
81  m.  Previous  tests  on  the  burner  had  indicated  X = 0.19,  and  data  from 
vertical  arrays  of  thermocouples  indicated  an  average  value  Ac  = 0.87.  Using 
all  of  the  above  data  the  test  run  is  found  to  be  characterized  by 
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m = 0,  e = 0.0032,  X = 1.3(109),  X£  = 

• . • 

Again,  assuming  the  wall  effect  to  be  significant  when  m^/m^  > 0.3  , i.e., 
(m^/m  )/X2  > 0.020  , figure  2 leads  to  the  conclusion  that  the  wall  effect 
must  be  taken  account  of  once  drops  to  0.85,  i.e.,  once  the  upper  layer 
thickness  exceeds  0.32  m. 

The  second  evaluation  of  the  test  series  will  involve  the  largest 

O 

constant  fire,  225  kW,  and  the  smallest  test  space  area,  40.6  m , of  the  test 

program.  This  configuration  had  a total  wall  perimeter  of  approximately 

43  m.  Tests  on  the  burner  indicated  Xr  = 0.24  for  burner  fires  in  the  range 

50  to  several  hundred  kW,  and  test  run  data  indicate  an  average  value 

X = 0.75.  All  of  this  leads  to 
c 

m = 0,  e = 0.027,  X = 1.1(1010),  X2  = 12. 

Using  the  previous  criterion  and  figure  2,  the  wall  effect  is  predicted  to  be 
significant  here  once  < 0.76. 

The  final  two  test  run  evaluations  will  refer  to  the  largest  and  smallest 
area  configurations  where  the  burner  was  controlled  to  generate  a fire  growing 
linearly  with  time  according  to 

Q(t)  = (t/2)  kW/s  (27) 

f 

Thus,  m = 1,  and  according  to  eq.  (20) 
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(28) 


Qc  = [ A/C2H3/2g1/2> ] kW/s 

First  the  largest  space  (A  = 89.6  P =*  81  m) : Here  Ar  = 0.24  , and 

test  data  indicate  an  average  value  A ^ = 0.81.  This  leads  to 

m = 1,  e = 0.0064,  X = 2.6(109),  X£  = 39. 

Using  the  previous  criterion  and  the  m = 1 .0  plots  of  figure  2,  the  wall 
effect  is  predicted  to  be  significant  once  < 0.93  . 

O 

Finally,  the  smallest  space  (A  = 40.6  m , P = 43  m) : Here  A = 0.24,  and 

test  data  indicate  an  average  value  A^  = 0.71.  This  leads  to 

tn  = 1,  e = 0.0063,  Xj  = 2.7(109),  X2  = 27. 

From  figure  2,  the  wall  effect  is  predicted  to  be  significant  once  < 0.90. 

5.  ESTIMATING  THE  WALL  TEMPERATURE 

Figure  2 and  the  above  example  calculations  were  based  on  the  assumption 
6 « 1 of  eq.  (26).  In  this  section  estimates  of  6 will  be  developed  so  that 

the  validity  of  this  assumption  can  be  evaluated,  and  so  that,  when  required, 

• • 

estimates  of  m^/m^  under  somewhat  more  general  nonambient  wall  temperature 
conditions  can  be  obtained  from  eqs.  (23)  and  (24). 

In  order  to  predict  6,  estimates  for  $ - 1 under  each  of  two  different 

wall  heating  conditions  will  be  established;  increases  of  <j>  - 1 due  to 
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1 at 


possible  radiation  from  the  fire's  combustion  zone,  and  increases  of  <f>w  - 

the  upper  wall  surfaces  due  to  convective  heating  from  the  wall  layer  flows. 

For  the  purpose  of  this  investigation  it  is  assumed  that  for  the  fire's  growth 

stage  of  present  interest,  and  far  enough  away  from  significant  ceiling  jet  - 

wall  flow  interaction,  one  of  these  two  mechanisms  will  dominate  the  increase 

of  d>  . 

T w 


5.1  Wall  Heating  Due  to  Radiation 


Increases  in  <J>W  due  to  radiation  may  be  the  more  significant  of  the 

heating  mechanisms  for  enclosure  spaces  which  are  configured  in  such  a manner 

that  a large  fraction  of  the  total  wall  surface  is  illuminated  by  the  fire's 

combustion  zone.  Under  such  circumstances  and  for  the  purpose  of  estimating 

<|>w  - 1 and  comparing  it  to  4>u  “ 1>  assume  the  fire  to  be  a characteristic 

distance  ( A/tt ) 1/2  from  all  wall  surfaces.  Then  the  average  rate  of  wall  heat 

transfer  due  to  radiation  from  the  combustion  zone,  q",  can  be  estimated  by 

R 


* Q(t) 

••  __  1 

qR  4A 


(29) 


For  the  purpose  of  the  present  estimate,  first  consider  wall  materials 
which  are  thermally  thick  up  to  some  time,  tj,  of  interest.  To  be  definite, 
it  is  reasonable  to  consider  a wall  material  as  thermally  thick  if  its  thick- 
ness, and  thermal  diffusivity,  <w,  satisfy 


LTHICK 


> 2/k  t_ 
w I 


(30) 
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For  tj-  of  3 minutes,  and  using  material  properties  from  table  1,  the  criterion 
would  be  satisfied,  for  example,  with  gypsum  board  walls  thicker  than  0.011  m 
or  concrete  walls  thicker  than  0.055  m. 


Using  the  thick  wall  assumption  and  eq.  (29)  in  Carslaw  and  Jaeger's 
solution  to  the  appropriate,  specified  flux,  heat  conduction  problem  for  the 
wall  leads  to  the  result 


w 


- 1 


X k1/2q 
r w o 

AT  Tfl/2Ak 

CO  w 


rt  t - n) 

j0  1/2 
n 


dn 


(31) 


where  kw  is  the  thermal  conductivity  of  the  wall  material. 


For  the  tm  fire  growth  studied  earlier,  eq.  (31)  can  be  integrated  to 

yield 


Radiation  to  Thermally  Thick  Wall: 


<(>  -1 
w 


1/4 

0.14X3Xj/ 


e(m+3)1/2r(m+l) 

T(m+l/2) 


(£ 

\0. 


_n/o  v 3(2m+l) 
^/j-l\2(m+3) 

140  / 


(32) 


6 = 


1/4 

0.43X3X^ 


(m+l)1/2r(m+l) 

r(m+l/2) 


(1-c. 


\ / 0.140  \ 


2(m+3) 


(33) 


where 


Xr  /l-Xr\5/4  (<\U2  /k  , 

“(1-V\1-V  VK/  \v 


Pr1/2H 

.1/2 


(34) 
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T(x)  i 8 the  Gamma  function,  Xj  is  defined  in  eq.  (25),  and  < and  k are  the 
thermal  diffusivity  and  conductivity  of  air  at  ambient  conditions. 

The  above  result  can  now  be  applied  to  the  previously  discussed  16.2  kW 
enclosure  fire  (gypsum  board  wall)  of  reference  2.  The  geometry  of  fire 
scenario  is  one  where  all  walls  of  the  enclosure  are  illuminated  by  the  fire, 
and  where  radiation  can  be  anticipated  to  dominate  the  wall  surface  heating. 
Using  the  previously  identified  parameters  of  that  fire,  and  the  k,  k values 
of  table  1 in  eq.  (33)  leads  to 

5 » (4>  - 1) / ( 4>  - 1)  = 0.29  at  C.  = 0.62 

w u l 

Note  that  this  ratio  is  small  enough  to  allow  the  6 ■ 0 estimate  of  figure  2 
to  be  useful.  Indeed,  if  the  approximation  of  eq.  (26)  is  supplanted  by  the 
newly  revised  estimate 

4>  - <{>  ° (<J>  ~ 1 ) ( 1 - 6)  » 0.71(<()  - 1) 

u w u u 

then  eqs.  (23)  and  (24)  would  predict  m^/m^  = 0.28  at  - 0.62  , instead  of 

• . • 

the  earlier  estimate  of  m /in  = 0.30  . 

w p 

If  the  wall  of  the  reference  2 enclosure  had  been  concrete  instead  of 
gypsum  board  (and  if  X would  not  have  been  significantly  different  from  the 
measured  value  of  0.84),  then  eq.  (33)  leads  to  the  estimate 

6 = (4>  - l)/(4>u  - 1)  = 0.09  at  C1  = 0.62 
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Besides  thermally  thick  walls,  wall  temperatures  for  thermally  thin  walls 
are  also  of  interest.  Again,  for  definiteness,  it  is  reasonable  to  consider  a 
wall  material  as  thermally  thin  if  its  thickness  satisfies 


lthin 


< 0.6/k  tT 
w I 


(35) 


Using  material  properties  from  reference  11  and  a t^  of  3 minutes,  this 
criterion  would  be  satisfied,  for  example,  with  mild  steel  walls  thinner  than 
0.028  m,  glass  walls  (windows)  thinner  than  0.006  m,  or  wood  panel  walls 
thinner  than  0.004  m. 

Using  the  radiant  heat  transfer  rate  of  eq.  (29)  and  the  tm  fire  growth 
leads  to 


Radiation  to  Thermally  Thin  Rear-Insulated  Wall: 


w 


TT  x4 


)3  ( m+ 1 ) 
(m+3) 


(36) 


6 = 0.25  x4(i  - e ) 


(37) 


where 


x4  -tt^t 


7 fc)  te)  (t) 


(38) 


If  the  wall  of  the  reference  2 enclosure  had  been  sheet  metal  construc- 
tion of  thickness  0.0015  m instead  of  gypsum  board  (and  if  X^  was  not  signifi- 
cantly different  than  the  measured  value  of  0.84),  then  eq . (37)  leads  to  the 
estimate 
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0.62 


6 » 0.13  at  » 

and  the  original  small  6 determination  of  the  significance  of  the  wall  effect 
continues  to  be  relevant. 


5.2  Wall  Heating  Due  to  Convection 

In  some  enclosure  fire  configurations,  radiant  heating  of  enclosure  wall 
surfaces  from  the  fire's  combustion  zone  may  play  a relatively  minor  role 
compared  to  convection.  This  would  be  the  case,  for  example,  if  direct,  line- 
of-sight,  combustion  zone  illumination  of  most  wall  surfaces  was  blocked  by 
furniture  or  segmented  partitions,  or  on  account  of  an  alcove-like  fire  loca- 
tion. Also,  in  a multi-room  enclosure  configuration,  which  is  an  extension  of 
the  alcove-like  single  room  configuration,  radiant  wall  heating  in  rooms  other 
than  the  room-of-f ire-origin  could  play  a minor  role  well  into  the  "smoke- 
filling” process. 

An  estimate  for  the  convection  driven  4>^( t ) increase,  to  be  obtained 
here,  will  be  based  on  an  eq.  (4)  heat  flux  to  the  wall  at  every  instant  of 

A tfc  ^ mmm 

time,  t , from  t = 0 to  t = t.  h(t  ) in  eq.  (4)  will  be  computed  from 
eq.  (5)  with  x = H - Z^(t*)  and  with  Grx  = Gr(t*)  according  to  eq.  (14)  or 
(24) . In  the  estimate,  5 will  be  neglected  compared  to  1 and  a laminar  wall 
layer  will  be  assumed  through  most  of  the  heating  history  up  to  the  time,  t, 
of  interest. 
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by 


All  of  the  above  considerations  lead  to  a convective  wall  flux,  q" 


— * 

q"(t  ) 


0.48  /m+3\5/4  

35/4H  VV  [l-?.]3/2 


15(m+l) 

4(m+3) 

T 


j|f 

where  t is  evaluated  from  eq.  (19)  with  t = t . Using  this  last  result 
reference  13  solution  to  the  specified  flux  heat  conduction  problem  for 
thermally  thick  wall  leads  to  the  result 


Convection  to  Thermally  Thick  Wall: 
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and  where  (x ;m) / 1^ (0 ;m)  has  been  evaluated  and  plotted  in  figure  3 for 

m = 0,  1.0,  2.0,  and  5.0. 


The  convective  wall  flux  of  eq . (39)  has  also  been  used  to  obtain 
$ - 1 and  6 for  thermally  thin  wall  constructions,  viz 


Convection  to  Thermally  Thin  Rear-Insulated  Wall: 
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and  where  I2(x ;m) / I2 ( 0 ;m)  has  been  evaluated  and  plotted  in  figure  3 for 
m = 0,  1.0,  2.0  and  5.0. 
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(49) 
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The  above  results  can  now  be  applied  to  the  previously  discussed  multi- 
room enclosure  fires  of  reference  1.  In  those  enclosure  configurations, 
radiant  heating  of  most  wall  surfaces  (i.e.,  outside  the  burn  room)  was  likely 
to  have  played  a minor  role  in  smoke  movement  phenomena  during  the  time  frames 
of  the  experiments.  Accordingly,  values  of  5 under  conditions  of  convective 
wall  layer  heating  have  been  computed  from  eqs.  (41)  and  (46).  These  are 
presented  in  table  2 along  with  the  results  of  all  previous  example  calcula- 
tions. As  can  be  seen,  for  the  cases  considered,  the  values  of  6 for  the 
convection  driven  wall  heating  examples  are  never  greater  than  0.011,  and  the 
original,  small  6 determination  of  the  significance  of  the  wall  effect 
continues  to  be  relevant. 

6.  QUASI STEAD I NESS  OF  THE  WALL  FLOWS  AND  A MEASURE  OF  THE 
LIKELIHOOD  OF  THEIR  PENETRATION 

The  validity  of  two  tacitly  assumed  aspects  of  the  wall  flow  will  be 
studied  in  this  section.  The  first  has  to  do  with  the  question  of  the  time 
required  to  reach  quasisteadiness  relative  to  an  assumed  steady  quiescent 
upper  layer.  Related  to  this,  the  second  has  to  do  with  the  expectation  that 
the  velocities  of  the  wall  flow  are  large  enough  to  actually  penetrate  the 
dropping  upper  layer  - lower  layer  interface. 

6.1  The  Quasisteady  Flow  Assumption 

The  time,  tg,  for  a wall  flow  to  establish  itself  as  a steady  flow 
(subsequent  to  a step  change  in  wall/ambient  temperature,  and  from  the  leading 
edge  of  the  boundary  layer  to  a station,  x)  has  been  studied  extensively  in 
the  literature.  Reference  7 presents  a review  of  this  work,  and,  in  terms  of 
the  present  nomenclature,  recommends  the  following  estimate  for  t 


25- 


ts  - 3(1  - Cl)2H2/(Gr1/2vJ 


(50) 


When  the  time  required  for  the  upper  layer  to  grow  to  a given  thickness 
and  temperature  is  at  least  as  large  as  the  above  t_  estimate,  then  there  is 

O 

good  reason  to  expect  that  the  quasisteady  wall  flow  assumption  is  valid. 

This  observation  together  with  the  previous  results  for  the  tm  enclosure  fire 
environment  finally  leads  to  the  following  criterion: 

With  regard  to  the  development  of  the  wall  flow  in  the  upper  layer  during 
tm  enclosure  fires,  the  assumption  of  quasisteadiness  is  justified  when  layer 
interface  elevations,  satisfy 


The  criterion  is  illustrated  in  the  upper  curve  of  Figure  4,  and  has  been 
successfully  used  to  test  the  validity  of  the  quasisteady  assumption  in  all 
previously  described  fire  scenarios. 


During  a fire  scenario  it  is  reasonable  to  expect  that  the  wall  flow  will 
penetrate  the  interface  when  the  characteristic  velocity  of  the  former,  say 
umax  ecl*  (10),  at  elevation  is  much  larger  than  that  of  the  latter, 
i.e.,  when 


(51) 


6.2  Wall  Layer  Penetration 


u, 


max 


(Z  = Z±)  » dZj/dt 


(52) 


-26- 


To  be  definite,  the  following  criterion  is  adopted: 


At  the  upper  layer  - lower  layer  interface,  penetration  of  the  wall  flow 
into  the  lower  layer  is  expected  if 

umax<Z  “ zl>  > 5dVdt 

Using  the  present  nomenclature  and  the  previous  results  for  the 
sure  fire  environment,  finally  leads  to  the  following  version  of  eq . 

/ -2/3  Al/2 

L_  ( _i ) > 1 Q HP_ 

?5/3  \ 0.140/  - yi-6J  yn+3 / AX2 

The  criterion  is  illustrated  in  the  lower  curve  of  Figure  4,  and  has 
successfully  used  to  test  the  validity  of  the  penetration  assumption 
previously  described  fire  scenarios. 

7.  SUMMARY  AND  CONCLUSIONS 

The  results  of  this  work  point  to  the  significance  of  buoyancy  driven 
wall  flows  on  the  development  of  hazardous  environments  within  enclosures 
containing  practical  growing  fires.  Such  flows  were  seen  to  be  important,  for 
example,  in  fire  scenarios  involving  relatively  weak  fires  and/or  plan  views 
with  relatively  long  peripheral  dimensions.  It  would  clearly  be  appropriate 
to  formulate  and  carry  out  an  experimental  program  for  quantitative  verifica- 
tion of  these  results.  In  any  event  it  appears  that  the  wall  flow  phenomenon 
must  be  accounted  for  in  mathematical  enclosure  fire  models  that  hope  to 


(53) 

tm  enclo- 
(53) 

(54) 

been 
in  all 
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estimate  dynamic  environments  generally  for  the  purpose  of  predicting  the 
response  of  fire  detectors  and  fire  intervention  hardware  and  the  time  of 


onset  of  conditions  detrimental  to  life  safety.  In  order  to  carry  out  such  an 
accounting,  further  analytic  and/or  experimental  research  on  various  aspects 
of  the  basic  buoyancy  driven  wall  flow  may  be  required,  and  should  be  pursued. 
However,  until  the  results  of  such  studies  are  available,  it  may  be  prudent  to 
add  wall  flow  algorithms  to  existing  models  on  a more  timely  basis.  Such 
algorithms  would  be  based  on  the  kinds  of  calculation,  estimates  and  results 
which  have  been  developed  in  the  present  work.  For  engineering  purposes,  such 
algorithms  may  prove  to  be  adequate,  indeed,  optimum,  even  in  the  long  term. 
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10.  NOMENCLATURE 


A 

area  of  enclosure 

Ci 

concentration  of  product  i 

S 

specific  heat  of  air 

Grx 

Grashoff  number,  eq . (2) 

Gr 

Grx  at  x = H - 

g 

acceleration  of  gravity 

H 

ceiling-to-f ire  distance 

h 

average  heat  transfer  coefficient 

*1*  l2 

integrals,  eqs.  (43)  and  (48) 

K 

kinematic  momentum  flux,  eq.  (8) 

k 

thermal  conductivity  of  air 

kw 

thermal  conductivity  of  wall  material 

lthick 

thermally  thick  wall  thickness 

lthin 

thermally  thin  wall  thickness 

m 

an  exponent 

• 

m 

mass  flux 

n 

an  exponent,  eqs.  (13)  and  (16) 

Nu 

X 

average  Nusselt  number,  eq.  (4) 

P 

perimeter 

Pr 

Prantl  number 

Q 

energy  release  rate  of  fire 

o * o 
O'  O' 

a characteristic  value  of  Q 
dimensionless  Qq,  eq.  (12) 

M 

q 

rate  of  heat  transfer  per  unit  area 

T 

absolute  temperature 

* 

t,  t 

time  from  ignition 
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a time  of  interest 


LI 

a Lime  ui  xutereuL 

u 

downward  velocity  in  wall  flow 

Xj  , ^2  , • • • 

dimensionless  parameters,  eqs.  (25),  (34),  (38),  (42), 
and  (47) 

X 

distance  below  ceiling 

y 

distance  from  wall 

z 

elevation  above  fire 

a 

a constant,  eq.  (16) 

r 

gamma  function 

6 

dimensionless  temperature  difference,  eq . (15) 

e 

dimensionless  parameter,  eq.  (18) 

C 

Z/H 

n 

dummy  variable 

K 

thermal  diffusivity  of  air 

K 

W 

thermal  diffusivity  of  wall  material 

X 

r 

fraction  of  Q lost  by  radiation 

X 

c 

fraction  of  Q transferred  to  enclosure  surface 

y 

viscosity 

V 

kinematic  viscosity,  y/p 

T 

dimensionless  time,  eq.  (19) 

<t> 

dimensionless  temperature,  T/T^ 

* 

Q/Qo 

Subscripts 


c 

convective 

i 

interface 

L 

lower  layer 

max 

maximum 
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p 

R 

u 

w 


plume 
radiant 
upper  layer 
wall  flow 
ambient 
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Figure  i.  Sketch  of  the  significant  features  ol  Llie  waLJ.  ellect 
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Figure  3.  Plots  of  (x;m)/ Ij (0 ; m)  and  I2(x;m)/ l2(0;m) 

per  eqs . (43)— (44)  and  (48)  — ( 49 ) , respectively,  for 
in  = 0,  1.0,  2.0  and  5.0 


= Zj/H 


1.0 

0.9 

0.8 


0.6 


Figure  4.  Plots  of  eq.  (51)  and  (54)  criteria  for  wall  layer 
quasisteadiness  and  penetration,  as  functions  of 
interface  elevation  and  for  tm  enclosure  fires 
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Table  1 . Thermal  Conductivity  and  Diffusivity  of 
Selected  Materials 


Material 

Conductivity 

Diffusivity 

m 

4-) 

air13 

0.024 

1 .9(10~5) 

gypsum  board^ 

0.17 

1.6(10“7) 

concrete  (1:2:4)^ 

0.92 

4.2(1 0—  7 ) 

steel  (0.1%  C)11 

46.0 

1.2(10-5) 
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Table  2.  Results  of  Calculations  on  the  Significance  of  the  Wall  Effect  in 
Full-Scale  Experiments  of  References  1 and  2 
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